On a heuristic point of view related to quantum nonequilibrium statistical mechanics 
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In this paper I propose a new way for counting the microstates of a system out of equilibrium. As, according 
to quantum mechanics, things happen as if a given particle can be found in more than one state at once, I extend 
this concept to propose the coherent access by a particle to the available states of a system. By coherent access I 
mean the possibility for the particle to act as if it is populating more than one microstate at once. This hypothesis 
has experimental implications, since the fhermodynamical probability and, as a consequence, the Bose-Einstein 
distribution as well as the argument of the Boltzmann factor is modified. 
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I. INTRODUCTION 



Nowadays, the development of new entropic forms has 
been followed by an increasing interest, as can be seen, for 
example, in Refs.|Ql Q]. Although it is possible to formulate 
new entropies from a strictly mathematical point of view [1], 
without connection with the physics implicated in such for- 
mulations, recently some works have appeared trying to un- 
derstand the link between the physical situation and the math- 
ematical formulation JH 0]. In this meantime, some works 
have appeared rising the question of a possible pseudononex- 
tensivity stemming from the generalized entropic forms J5i0]. 

In general, nonextensive formulations are related to 
nonequilibrium situations, where the Boltzmann factor 
exp(— @E), presumably, plays not a preponderant role, being 
difficult, if not impossible, to associate a definite temperature 
to the system. In some cases, however, by considering sit- 
uations only slightly out of the equilibrium, it is possible to 
ascribe a temperature to the system, which results in a distri- 
bution function different from that of Boltzmann J7[]. 

In this paper, inspired by some ideas from quantum optics, 
I propose a new way for counting accessible states to a given 
particle, in such a way that its thermodynamical probability 

is modified, with direct consequences in the entropic form 
S oc £1 of the system. As is well known in the quantum op- 
tics domain, which deals fundamentally with nonequilibrium 
systems, an initially pure state can be described, in its most 
general form, as a superposition of each state physically ac- 
cessible to the particle. The role of the reservoir, even at the 
idealized zero temperature, is to lead the system to a com- 
plete mixture at the end of the so-called decoherence time td- 
Thus, even before the thermalization occurs the loss of coher- 
ence of the system, or, in other words, the system capacity to 
access, coherently, every possible state. This state of affairs 
suggests an entirely new way to count the accessible states. It 
is this connection, until now not explored, between the new 
way to count the accessible states and its consequences to the 
the entropic form of the system that we will explore in the next 
sections. 



II. REDEFINING THE MICROSTATES 

In statistical mechanics, to define a microstate it is neces- 
sary to take into account the (un)distinguishability of the par- 
ticles, which gives rise to different configurations (see Tab. 1). 
For calculating all the possible configurations we now take 
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TABLE I: The configuration of two accessible states for a) two in- 
distinguishable particles and b) two distinguishable particles. (1) de- 
notes the first available state and (2) denotes the second available 
state. 

into account, beside this characteristic, this another one: the 
possibility to the particle simultaneously access more than one 
state, or, to avoid eventual difficulties related to interpretations 
matter inherent to the quantum formalism, the possibility to 
the particle to coherently access the available states. This sit- 
uation is shown in Tab. 2 for the case of two identical particles 
having two accessible states. Note that if the particles are dis- 
tinguishable, the corresponding configuration is different. 
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TABLE II: A system out of equilibrium composed by two particles 
having two accessible states. (1) denotes the first available state, (2) 
denotes the second available state, and (12) denotes the coherent ac- 
cess to both states. 
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Comparing Tab.I and Tab. II, we see that, clearly, the 
nonequilibrium situation requires a new way for counting mi- 
crostates. This new way to count, shown in Tab.2, can be rep- 
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resented by the following sequences, where the number be- 
tween parentheses indicates the state occupied and the letter 
following the parenthesis indicates the corresponding occupa- 
tion by the particle a, which is identical to all the others: 



(l)o(2)o(12); (l)o(2)(12)a; (l)(2)o(12)o; 
(l)aa(2)(12); (l)(2)aa(12); (l)(2)(12)aa. 



(1) 



As for example, the last sequence, (l)(2)(12)aa, corresponds 
to two particles accessing coherently the two states (1) and 
(2), while the second sequence, (l)a(2)(12)a corresponds to 
one particle accessing the state (1) and the other accessing 
coherently the states (1) and (2). Note that, as the sequence 
must initiate by a number, and existing three possible number 
of states, 1,2, and 12, will remain 3—1 numbers plus two 
letters a (particles) to be set in whatever order (permutation). 
Therefore, the number of unrepeated sequences is 



3 x (3 - 1 + 2)! 
2!3! 



= 6, 



(2) 



where we have put a superscript (*) to remind us that we are 
treating with nonequilibrium situation. Proceeding in a gen- 
eral manner, for gj sublevels with N* particles, the number 
w* of unrepeated sequences is 



Gj(Gj + N* - 1)! _ (Gj + N* - 1)! 

~ (Gj-1)\N*\ ' 



GjWp. 



(3) 



where Gj = Ysk=i Ggj.k * s tne num ber of possible sequences 
formed from gj, and C n ^ m — n\/(n — m)\m\. Taking as 
example the configuration given by Tab.2, where gj = 1, 



then 



2, Gj = ELi c 2,fc, thus Gj = C 2 



c 2 



(3 
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6, 



3; 



(4) 



which is the number of sequences given in Eq.([T]i correspond- 
ing to Tab.2. Therefore, the nonequilibrium thermodynamical 
probability wl for a given macrostate k is 



n 



(Gj + N*-1)\ 



(Gj 



i)\m\ 



(5) 



As G. 



_ \^9j 



i G gj ,k 



= C, 



Efe= 2 C g ./.aiuir, =g jf 



letting L aj =Ek=2G gj ,k 



then Eq.© can be written as 



n 



( gj + L gj + N;-iy. 

(9j 



L gj - l)\N*\ 



(6) 



From Eq.© we can see that the only changing in the ther- 
modynamical probability is the appearance of the factor L g j 
modifying the degeneracy gj, and, as a consequence, modi- 
fying also the number of macrostates, £1 = J^k w k> an< ^ tne 
entropy of the system. Before ending this section, we call 
attention to the plausibility in presume that, given a system 
out of the equilibrium with N* particles and n levels, each 



of this having gj sublevels, as the equilibrium is established 
(L g j — > 0), the N* particles of the system accommodate by 
the rij levels, with each level receiving N* particles, which 
are distributed by the sublevels. Also, as it is easily verified, 
Eq.© gives rise to a Bose-Einstein-like statistics, with gj re- 
placed by Gj . That this is so can be checked in the follow- 
ing manner, proceeding by analogy with the equilibrium situ- 
ation: First, we take the In from both sides of Eq.©. Second, 
we use the Stirling formula. Third, we differentiate with re- 
spect to Nj and use d In w*/ dN* — e* , where e* generalizes 
Ej = f3Ej, (3=1 /kT , to find 



N*/Gj 



exp(e*) - 1 



(7) 



This is an interesting point, having experimental implication: 
the Bose-Einstein statistics is corrected, since the equality 
Gj = gj and e* = £j will be valid only when the complete 
equilibrium is reestablished. Thus, for systems only slightly 
out of the equilibrium, the energy emitted should be slightly 
different from that corresponding to the system in equilibrium. 
Note that, as e* = Ej = 0Ej when the equilibrium is restated, 
it is convenient to expand e* in power series of e 



de* 



2! de 2 , j 



1 dh* 3 
3! de* £j ' 



(8) 



which, requiring that e* — > ej = j3Ej when the equilibrium 

is restated, gives £q = and ^ = 1, such that the first order 
correction to the Bose-Einstein distribution can be explicitly 



written as 



N*/Gj 



1 



exp PE + ai ((3E)' 



2! de 2 t 



Oil. 



where we have kept only a few terms and put 

Note that from this approach the net effect stemming from 
the nonequilibrium on a given system is the increasing in 
the degeneracy, which in turn increases the available states 
given by Q. The Boltzmann factor, to be recovered when 
exp (/?*e*) ^> 1 , is modified, and we will explore more about 
this in the next Section. The choice of the more convenient 
entropic form associated to this new thermodynamical proba- 
bility is discussed in the last Section. 



III. THE NONEQUILIBRIUM PARTITION FUNCTION 

Let us focus our attention to the bosonic particles, since the 
other cases are similar. By definition, the partition function is 
defined as a sum in all microstates (ms): 



Z = £)exp(-/3E), 



(9) 



where E is the energy of the system and j3 is related to the 
temperature T of the system by the Boltzmann constant [3 
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= 1/ kT. Writing the energy E in terms of the number of 
particles in the state of z — th energy e(i) of the system, we 

will have E = ^ i n i e(i). Of course, in this case the total 
number of particles is simply N = J^i n i- 

For an out of equilibrium system, we introduce the co- 
herent access hypothesis to several states, which consists in 
maintaining the same form as that of Eq.©, but replacing 
J2i rne(i) by V, ; mje(i,j, ...), where must be inter- 
preted as being the number of particles coherently accessing 
the energy levels e(i) and e(j). For example, as discussed in 
Section I and represented in Tab.2, e(i,j) represents the co- 
herent access related to the energy levels i and j, and e(l, 2) 
represents, for example, the states (1) and (2) being coher- 
ently populated. 

For demonstrating that the partition function preserves its 
form given by Eq.© even at the nonequilibrium situation, it is 
enough to maintain this following postulate, which is valid for 
equilibrium situation: that two systems, in contact with a third 
one, as for example a reservoir at temperature T, act indepen- 
dently of each other while both the systems exchange energy 
with the reservoir. Although this demonstration is straightfor- 
ward, for completeness we address the reader to the appendix. 
Continuing to denote the nonequilibrium quantities with a su- 
perscript (*), thus according to Eq.d26b of the appendix, if 
P(e* = (3*E*) is the probability for a given system out of 
the equilibrium is in a particular microstate whose configura- 
tion is described by e* = (3*E*, then 



P(e*) = 



exp(-e*) 



(10) 



Now, using Eq.® and requiring that e* — » (3Ej when the 
equilibrium is restated, the EqdTOb can now be written as 



P(e*) = — exp -PEj-aitfEj) 



a 2 {pEjf + ota CU\,, 



(ID 



where the other constants were renamed for convenience as 
^j -gpr- = oc n -\. Such a state of affairs giving origin to an 
infinite number of free parameters was studied in Refs. 
in a different context. Note that for systems only slightly out 
of the equilibrium this last equation can be written as 



P{E) = — exp \-/3E - ati (f3E) 2 



(12) 



where we have dropped out the superscript and the index i. 
Some experiments seem to point for the importance of this 
last term, which modifies the Boltzmann factor Q. 



IV. CONNECTION WITH ENTROPIC FORMS 

As discussed in Section I, since the thermodynamical prob- 
ability was modified, a natural question emerging is what is 
the best entropic form related to it. Of course, depending 



on our choice we will face with different implications. Once 
there is a plenty of entropic forms at our disposal, we will fo- 
cus our attention only in two of them: the Boltzmann-Gibbs 
(Sbg) an d the Tsallis (S q ) entropies. As is well known, while 
the first is extensive, /.e. S B g{A + B) = S B g{A) + Sbg(B), 
the second in general is not, i.e., S q (A+ B) ^ S q (A) + S q (B) 
if'/ * I- 

Let us begin adopting the Boltzmann-Gibbs entropy, as- 
suming for now that the single effect of the nonequilibrium 
is to increase the degeneracy of the system, as seen in Sec- 
tion I. It will be possible to reconcile Eq.dTTb to an extensive 
entropic form such as that of Boltzmann and Gibbs? Indeed, 
that this is possible was shown in Ref.[5], in the following 
way. Given the density operator p of the system and the Boltz- 
mann constant k, for maximizing the Boltzmann-Gibbs en- 
tropy Sbg = —kTrplnp subjected to the constraints given 
by the moments 



{(AE) n ) = TrpH n , 



(13) 



n integer, we vary p in Sbg an d in those for the constraints, 
Eq.([T3l, multiplying each constraint by the undetermined La- 
grange multiplier f3 n , and adding the result, obtaining 



TV ( 1 + PnH n + lnp Sp = 0. 



(14) 



Since all the variations are independent and dp is arbitrary, 
it follows the extended (non-Maxwellian) distribution In p = 

oc 

— 1 — PnH n , or, equivalently 

n=0 



oo 



p= Z- L cxp(-} B n H n ) 



(15) 



where the partition function is Z = Trexp(— ^ (3 n H n )- In 

71=1 

the energy representation where H\E) = E \E), Eq.(fL5b now 
reads, 

oo 

P(E) = Z- 1 exp(- PnE n ) = Z- 1 x 

71=1 

exp (-faE + (3 2 E 2 + (3 3 E 3 + (3 4 E\..) (16) 

oo 

with Z = X)_E ex P( — S PnE n ). The Lagrange multipliers 

71=1 

Pk 316 formally obtained from /3k = — g^if , considering E k 
— Yfe as independent variables. The equality between Eq.(fT6b 
and Eq.lfTTll is guaranteed, provided that n = a n _i/3™ and 
(3± = ao = (3. Therefore, according to this view nonequilib- 
rium systems remains extensive, although requiring a poste- 
riori knowledge of the variance (second central moment), the 
coefficient of skewness (third central moment), the kurtosis 
(fourth central moment), and so on, thus giving rise virtually 
to an infinite number of free parameters. 

Of course, instead of using infinite parameters, we could 
just use a single one by redefining a new ensemble fully de- 
termined by this single parameter. An aesthetically appealing 
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way to do so is to expand Eq.dTTb in terms of the Tsallis en- 
tropic index [8], as we will see in a moment. Consider the 
following expanded form of Eq. dTTb : 



cxp —ftEj — 



my 



(17) 



where in general a„ 



(<?-i) r 



This is equivalent to the 



statement that the old ensemble which depended of 0, {a n } 
and Ej becomes now a function of only j3, q and Ej, Eq.(fT7t 
can be rewritten as 



1 



P(Ej) = -exp 



1 



(l-q)(3E 3 



[fiEiY- 



m 3 y 



(18) 



where it is easily recognized the expanded form of the log- 
arithm function ln(l 



x = (1 — q)/3Ej, such that Eq.([T8l becomes 



P(E) = -[!-(! -q)f3E 3 ] — 



(19) 



which is the g-distribution stemming from the extremization 
of Tsallis entropy 



1 



Sg — k- 



(20) 



when considering a family of constraints determined by the 
g-expectation value of the energy 



(21) 



besides the norm constraint YlPj = 1- Therefore, a formal 

j 

agreement between Tsallis and Boltzmann-Gibbs entropies is 
afforded. As pointed out in Ref.||5|], this formal equivalence 
between the Boltzmann-Gibbs and Tsallis entropy gives rise 
to an important issue related to a possible pseudononextensiv- 
ity. 



consequence of the coherent access hypotheses, the process 
of counting the possible states of a physical system is modi- 
fied. I have found a modification on both Bose-Einstein and 
Boltzmann-Gibbs distribution, which is in principle experi- 
mentally detectable. Actually, it is possible that the correction 
to the Boltzmann factor obtained by the method developed 
here is the one suggested by some experiments [7]. Although 
I have exemplified for the specific case of bosons, the exten- 
sion to fermions is straightforward. Finally, I expect that the 
coherent access hypothesis introduced here eventually makes 
possible the exploration of new ways of treating problems re- 
lated to nonequilibrium situations, or differing from the equi- 
librium in a slightly manner. 



Appendix I 

To demonstrate that the partition function and the Boltz- 
mann factor retain the same form as Eq.© in the nonequi- 
librium situation, it is enough to follow the usual derivation, 
as for example, that given in Ref. lfToll . Thus, consider a sys- 
tem composed by two subsystems A and B. The probability 
for this composed system to be in the energy state E* A+B is 
P A + B (E A+B ), where the superscript (*) remind us that the 
system is out of equilibrium. If, as usual, the interaction en- 
ergy can be neglected, thus the energy of the composed system 



is E A+B = E *A 



E* B , and 



P A+B {E* A+B )=P A {E* A )+P B {E* t 



(22) 



is the probability for the composed system to be in a particu- 
lar state such that the subsystem A has an energy E A , and, at 
the same time, the subsystem B has an energy E* B . Now, sup- 
pose that these two subsystems is put in contact with a third 
system, for example, a reservoir at temperature T. While per- 
sisting the nonequilibrium situation (and even after that), the 
two subsystems A and B act independently of each other, with 
both subsystems eventually exchanging energy with the reser- 
voir. Beside that, the energy exchanged with the reservoir by 
a given subsystem does not influence the energy that the other 
subsystem can exchange with this same reservoir. This as- 
sumption, valid for two systems in equilibrium with a reser- 
voir, is here assumed to be valid also when the equilibrium 
was not reached. Therefore, as these events are independent, 
we can write 



P{E A+B ) = P(E A )P(E* B ). 



(23) 



Differentiating Eq.(|23l with respect to E* A and E* B and equat- 
ing this result we obtain (dP/ dE* = P') 



V. CONCLUSION 

In this paper I explored an analogy between the nonequi- 
librium thermodynamics and some well-established situations 
from quantum optics, concerning the problem of coherent ac- 
cess to the multiple states available to a given particle. As a 



P' A {E A )P B {E B ) = P A (E\)P' B {E B ). (24) 

Next, separating the variables and equating the result to a con- 
stant, we have 



P'a(E a ) _ P' B {E B ) 



P A (E* A ) P B (E* B ) 



(25) 
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where (3* is a constant independent from either E* A or E* B . Of 
course, in the equilibrium situation we must have = 
1 / kT. From Eg. (1251) follows, therefore, our desired result 

, „ s exp(-f3*E*) 
P(E*) = — ^ -, (26) 

where the partition function for the nonequilibrium situation 
is Z* and the index were dropped given the validity of Eq.(l26t 



for the two subsystems. 
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